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1. Introduction 

A deeper understanding of the mathematical foundation of phenomeno¬ 
logical models of the atomic nucleus is often related to the studies in simple 
models for which rigorous solutions can be obtained. In the domain of 
weakly bound or unbound exotic nuclei, we are missing such models, even 
though it is well known that the coupling to continuum may enhance the 
pairing correlations and the stability of weakly-bound nuclei [1]. 

The objective of this Letter is to derive an approximate solution for the 
schematic pairing Hamiltonian in the single-particle (s.p.) space involving 
both bound and continuum states. This model is a generalization of the 
rational Richardson-Gaudin (RG) model [21 [3] in the Berggren ensemble |3] 
of s.p. states. We shall generalize the solution found by Richardson for the 
rational RG model to describe the pairing interaction of particles in bound 
and continuum s.p. states. 

2. Pairing Hamiltonian in the Berggren basis 

The pairing Hamiltonian for the rational RG model is given by: 

D Q D 

H = Y^ eaclca + jYl cl4c,3C/3 (1) 

a 0,0 
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where ea and G are the energies of (bound) s.p. states and the pairing 
strength, respectively. Here, D stands for the number of s.p states, i.e 
D = 2N + V with N the number of pairs, u = the number of unpaired 

particles, and U the number of levels. The operators ct (cq) stand for the 
particle creation (annihilation) operators and a = {a, uIq,} = {na,ia,ja,rna}- 
It is convenient to rewrite this Hamiltonian in terms of pair creation 
(annihilation) operators b]^{ba) and particle number operator ha ■ 

= , bl = 1/2 Y, 44 = (fta)^ (2) 

ma. ma 

Then, the Hamiltonian becomes : 

u u 

H = Y X! 4^a' , (3) 

a,a' 


Expressing the pairing Hamiltonian ([T]) in the complete one-body Berggren 
basis [3], including bound (6), resonance (r), and non-resonant (c) states, 
one obtains: 



Note that the Hamiltonian of Eq.(I3|) is not exactly the same as in Eq.([T]), as 
the latter is defined in the Hilbert space, whereas the former is defined in 
the rigged Hilbert space, comprising bound, resonant, and scattering states 
[5]. In this equation, the summation over c {c') represents the summation 
over different partial waves {£,j) until (^max, jmax)- is related to the s.p 
energy: Cc = h?ka/2m, with m the particle mass. The discrete sums run 
over bound s.p. states and s.p. resonances inside the complex plane between 
the fc-contour L~^ and the real fe-axis. More about the complete Berggren 
s.p. basis and its application can be found in Ref. [B]. The operators 
satisfy the commutation relations: 


ni, 4 

bi,b\, 


‘25iiibl 



(5) 
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for discrete (bound and resonance states) s.p states, and 




26 {h - k',)bl 

6{kc - k'c)—^ ± bk„k'„fik„ 


( 6 ) 


for the non-resonant scattering s.p states. Here, ^i^k„ is the degeneracy 
defined as ^i^kc — ‘^ji,kc + 1- A Kronecker delta involving scattering states 
had to be introduced in Eq.Q to preserve antisymmetry. The upper (lower) 
signs in these equations stand for bosons (fermions). 

To discretize continuum in Eq. 0, it is convenient to define new oper¬ 
ators: b}j = b]jyfw^ and fiq = Wqhq , where Wq is a weight which equals one 
for bound and resonance states, and is equal to the Gauss-Legendre quadra¬ 
ture weight for scattering states. Index q runs over bound, resonance and 
discretized scattering states which are all normalized to one. New operators 
satisfy: 


n, 






qq "q 


The pairing Hamiltonian expressed in these operators: 


u 


u 


H — ^2 ^q^q + X! ^qq'^q'^g ’ ^qq' ~ 


q,q 


contains the state-dependent pairing force and is not integrable [7j. 


(7) 

( 8 ) 

(9) 


3. The approximate solution of the pairing Hamiltonian 

An approximate solution for the pairing Hamiltonian ([9]) in the contin¬ 
uum can be found by replacing the Kronecker delta in the commutation 
relation ([U]) by the Dirac delta : 




26 {kc 



With this change, the new pair operators 6j(6q) satisfy: 


bq^K, 


= 26, 


qq 



nq 

2Wq 


( 10 ) 


( 11 ) 
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and the commutation relation involving remains the same. 

Let us now derive the solution of a Schrodinger equation for fermions: 


N 


U 


L^l^n:norm >— 


n ^ ninorm 


> ; En — ^ Ej^ + ^ Vne, 


q^q 


( 12 ) 


7=1 


where Ej^ are the pair energies. The derivation for bosons is similar and 
will not be given here. 

For the many-body state, we take the following ansatz : 


IT 


N 

n;norm) = -®,7.^;norml^) > '®J-,;norm 

7=1 


U 




, “1 


•^q V ^9 

2eg — Ej 


where the normalization constants cj are given by: 


u 


{cjM) Cj^ 


= E 


Wn 


{2eq - Ej^ 


Here |v) is a state of unpaired particles satisfying: 

bq\l^)=0 , hq\u) = iyq\v) 

To simplify, we define and write: 


IT 


N N 

n;norm) = H =Cn\^n) ; Cn = H Cj^ 

7=1 7=1 


Let us begin by evaluating the commutator: 


N 


n < 


7=1 


N I /7-I 


N 


=E n^i n 


7=1 \ \? 7=1 


VM=7+1 


For that, we rewrite the Hamiltonian (j3|) in a discretized form: 

u u 

■^ = E ^9^-? + GbIBq ; ^ b\y/EEi , 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


where U is the total number of bound, resonance and discretized continuum 
states. Knowing the commutation relations for bq,bq, and hg, one obtains: 


t 


hg,Hy 


2eg - Ej 


Bo,B\ 


U 


= E 

9 


Wq^lqj2 - flq 


2eg - Ej^ 


(19) 
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and 


N 


n, 


n B] 

^=7+1 


N 


= E 


/j'=7+l 


{2eq - Ej ) \ 


N 


n 


m=7+1;^m' 


( 20 ) 


Then, using Eq.([T9|) we get : 


H,B 


t 




;t ^ r1 /1 I T 


= Ej,B'j +BUI+GY, 


2 eg — Ej^ 


( 21 ) 


Inserting (12X1) and (l2n]i into Eq. m, one finds: 




B, n M 


7=1 


B, 


t 


N 


= E^tI -^EE2e -Ej 

7=1 q \7=1;^7 


n 


Un 


7=1 

N / U 

+ E 1 + E 

7=1 \ 

N U 

- EE 


7=1 


N 


GWi^q/2 ^ I Y[ b] 


2eq — Ej 


q 


, 7=l;^7 


2G^Blbl 


N 


n 


p=J+l 


( 22 ) 


Now, applying the commutator (|22p on the vacuum state \v), and using the 
relation: 

u 


E 


-bI 


{2eq - Ej^){2eq - Ej^) Ej^ - Ej^ 


(23) 


we obtain: 


«|W„,. .I:' (:. E ( n' 4. 


-EE 


7=1 \ q 

2G 


^7=l;^7 


7=1 \/x=l;7^7 Bj^ 


^pair 

41 n 4, 

I7=l;^7 


(24) 


As H\^n) = En\^n)^ One obtains the generalized Richardson equation for 
pair energies of the discretized pairing Hamiltonian in the Berggren basis: 


u 


N 


_ GWq {Vq Hg/2) 


2G 


Bji Bj^ 


= 0 


2eq — Ej^ 


(25) 
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In the continuum limit, when the number of discretized scattering states 
becomes infinite and their weights go to zero, these equations take the form: 


1-E 


2Gdj 


ieb,r 


^max ijmax 

- E 


2.Gdk^ 


N 


'l+ h^k^/m — Ej. 


-dkc + ^ 


2G 




= 0 


with di = Vil2 — Vli/A. 


(26) 


4. Discussion and conclusions 

It is interesting to discuss limiting cases of the generalized Richardson 
equation ([25]). If we deal with a discrete set of bound s.p. levels, the weights 
Wq are all equal one and Eq. ([2^ reduces to the exact Richardson solution 
for the rational RG model mm- By the same argument, Eq. (|25l) provides 
an exact solution of the pairing model with a continuum treated in the pole 
approximation, i.e. neglecting the non-resonant continuum states. Einally, 
Eq. (j25D is also the exact solution of the pairing model in the non-resonant 
continuum because in this case one may take the same weights Wq = w 
for all states q and renormalize the pairing strength G' = Gw. At this 
level we no longer use the Gauss-Legendre quadrature, as we start from the 
Richardson equations taken at continuum limit, so that it is possible to use 
another discretization scheme where all weights are equal for simplicity. In 
this particular case, the third sum in Eq. ([25]) goes to zero and one obtains: 


G 


dk 


2efc — Ei 


-dk = 1 . 


(27) 


In a most general case of the rational RG Hamiltonian in the Berggren 
ensemble of s.p. states, Eq. (I2^ provides a reliable approximation of the 
exact solution in the limit of weak pairing correlations. The quantitative 
comparison between results of Eq. (j25p and exact results of Gamow Shell 
Model mm for the pairing Hamiltonian in the Berggren basis will be 
discussed elsewhere m- 

We wish to thank J. Dukelsky for useful discussions. This work was 
supported in part by the GOPIN and GOPIGAL French-Polish scientific 
exchange programs. 
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